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Abstract. In this paper our aim is to show some mean value inequalities for the modified Bessel 
functions of the first and second kinds. Our proofs are based on some bounds for the logarithmic 
derivatives of these functions, which are in fact equivalent to the corresponding Turan type inequalities 
for these functions. As an application of the results concerning the modified Bessel function of the 
second kind we prove that the cumulative distribution function of the gamma-gamma distribution is 
log-concave. At the end of this paper several open problems are posed, which may be of interest for 
*w ^ further research. 
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Let us consider the probability density function ip : M — > (0, oo) and the reliability (or survival) function 
$ : R — > (0, 1) of the standard normal distribution, defined by 

1 -i /*oo 

p{u) = -^e- u ' 2 / 2 and *(«) = -=/ e~ t2/2 dt. 

V27T V27T J u 

The function r : R — > (0, oo), defined by 



< 

u 

^: r(tl) = !M = e ^r e ^ dt , 

Ju 

is known in literature as Mills' ratio [23 sect. 2.26] of the standard normal distribution, while its reciprocal 
1/r, defined by \/r(u) = <p(u)/$(u), is the so-called failure (hazard) rate, which arises frequently in 
economics and engineering sciences. Recently, among other things, Baricz |12[ Corollary 2.6] by using 
the Pinelis' version of the monotone form of l'Hospital's rule (see [3SJ 121 H] for further details) proved the 
. following result concerning the Mills ratio of the standard normal distribution: 



Theorem A. If u\,u 2 > Uo, where uq rs 1.161527889 . . . is the unique positive root of the transcendent 
equation u(u 2 + 2)$(u) = (u 2 + l)p(u), then the following chain of inequalities holds 

OV / -i \ 2r(m)r(u 2 ) ( Ux + u 2 \ n — r-p — r , r( Ul ) + r(u 2 ) ( 2 Ul u 2 

Moreover, the first, second, third and fifth inequalities hold for all Ui,u 2 positive real numbers, while the 
fourth inequality is reversed if u\,u 2 € (0,uq). In each of the above inequalities equality holds if and only 
ifui = u 2 . 

X 

We note here that, since Mills' ratio r is continuous, the second and third inequalities in (QJ) mean 
actually that under the aforementioned assumptions Mills' ratio is log-convex and geometrically concave 
on the corresponding interval. More precisely, by definition a function / : [a, b] C R — > (0, oo) is log-convex 
if In / is convex, i.e. if for all ui,u 2 € [a, b] and A £ [0, 1] we have 

/(A«i + (i - \)u 2 ) < [/Ml'l/M] 1 ^. 

Similarly, a function g : [a, b] C (0, oo) — > (0, oo) is said to be geometrically (or multiplicatively) convex 
if g is convex with respect to the geometric mean, i.e. if for all u\,u 2 € [a, b] and A G [0, 1] we have 

jK^^bK^fe)] 1 ^. 

We note that if / and g are differentiable then / is log-convex if and only if u i— > f'(u)/ f(u) is increasing 
on [a, b], while g is geometrically convex if and only if u i— > ug'(u)/g(u) is increasing on [a, b]. A similar 
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definition and characterization of differentiable log-concave and geometrically concave functions also 
holds. 

Mean value inequalities similar to those presented above appear also in the recent literature explicitly 
or implicitly for other special functions, like the Euler gamma function and its logarithmic derivative 
(see for example the paper [2] and the references therein), the Gaussian and Kummer hypergeometric 
functions, generalized Bessel functions of the first kind, general power series (see the papers [SHHH2], and 
the references therein), Bessel and modified Bessel functions of the first kind (see [13l H8l 132] ). 

In this paper, motivated by the above results, we are mainly interested in mean value functional 
inequalities concerning modified Bessel functions of the first and second kinds. The detailed content 
is as follows: in section 2 we present some preliminary results concerning some tight lower and upper 
bounds for the logarithmic derivative of the modified Bessel functions of the first and second kinds. These 
results will be applied in the sequel to obtain some interesting chain of inequalities for modified Bessel 
functions of the first and second kinds analogous to ([1]). To achieve our goal in section 2 we present 
some monotonicity properties of some functions which involve the modified Bessel functions of the first 
and second kinds. Section 3 is devoted to the study of the convexity with respect to Holder (or power) 
means of modified Bessel functions of the first and second kinds. The results stated here complete and 
extend the results from section 2. As an application of our results stated in section 2, in section 4 we 
show that the cumulative distribution function of the three parameter gamma-gamma distribution is 
log-concave for arbitrary shape parameters. This result may be useful in problems of information theory 
and communications. Finally, in section 5 we present some interesting open problems, which may be of 
interest for further research. 

2. Monotonicity properties of some functions involving modified Bessel functions 

As usual, in what follows let us denote by l v and K v the modified Bessel functions of the first and 
second kinds of real order v (see 41 J, which are in fact the linearly independent particular solutions of 
the second order modified Bessel homogeneous linear differential equation [JTJ p. 77] 

(2) u 2 v"(u) + uv'(u) - (u 2 + v 2 )v(u) = 0. 

Recall that the modified Bessel function l v of the first kind has the series representation [411 P- 77] 



{u/2) 2n +» 
n\r(n + is+ 1)' 



n>0 

where v ^ —1,-2,... and u € K, while the modified Bessel function of the second kind K u (called 
sometimes as the MacDonald or Hankel function), is usually defined also as 41, p. 78] 

= 

2 sin vtt 

where the right-hand side of this equation is replaced by its limiting value if v is an integer or zero. We 
note that for all v natural and a€l we have I v (u) = I— U (u), and from the above series representation 
I v {u) > for all v > — 1 and u > 0. Similarly, by using the familiar integral representation [JTJ p. 181] 

/>oo 

(3) K„{u)= e - ucosh *cosh(z4)dt, 

Jo 

which holds for each u > and fEl, one can see easily that K v (u) > for all u > and 

The following results provide some tight lower and upper bounds for the logarithmic derivatives of 
the modified Bessel functions of the first and second kinds I u and K v and will be used frequently in the 
sequel. 

Lemma B. For all u > and v > the following inequalities hold 

(4) I^ U 2 +V 2 < ^IM. < ^2^2". 

w V v+ 1 I v (u) 

Moreover, the right-hand side of Q holds true for all v > — 1. 



Lemma C. For all u > and V > 1 the following inequalities hold 
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Moreover, the right-hand side of ([5]) holds true for all v 6 



The left-hand side of (j4]) was proved for u > and positive integer v by Phillips and Malin [34], and 
later by Baricz [14] for u > and v > real. The right-hand side of (j4]) appeared first in Gronwall's 
paper [27, for u > and ^ > (motivated by a problem in wave mechanics), it was proved also by 
Phillips and Malin [33] for u > and v > 1 integer, and recently by Baricz [T3J for u > and v > —1/2 
real (motivated by a problem in biophysics; see [53]). For this inequality the case u > and ^ > — 1 real 
has been proved recently in [T7] . 

The left-hand side of ([5]) was proved first by Phillips and Malin 34 for u > and v > 1 positive 
integer, and was extended to the case u > and i/ > 1 real recently by Baricz [T7]. Finally, the right-hand 
side of ([5]) was proved first by Phillips and Malin [33] for u > and i> > 1 integer, and later extended to 
the case of u > and v real arbitrary by Baricz [T3] . 

It is worth mentioning that the inequalities (jU) and ([5]), which have been proved recently also by Segura 
[37], are in fact equivalent to the Turan type inequalities for the modified Bessel functions of the first 
and second kinds. For further details the interested reader is referred to [131 H7J HU [3D] [37] and to the 
references therein. 

Our first main result reads as follows. 

Theorem 1. The following assertions are true: 

(a) u i — y ul' v (u) I I„(u) is strictly decreasing on (0, oo) for all v > 1; 

(b) u I—?- ul'yiu) / I v {u) is strictly increasing on (0, oo) for all v > —1; 

(c) u i— > y/ul v (u) is strictly log-concave on (0,oo) for all v > 1/2; 

(d) u i — ^ u 2 'I' v (u) I ' Iy{u) is strictly decreasing on (0, oo) for all v > i/q, where 
the positive root of the cubic equation 8^ 3 — %v 2 — 2v — 1 = 0. 

In particular, for all Ui,U2 > and v > v q the following chain of inequalities holds 



1.373318506. 



(G) 



2I„(ui)I v {u 2 ) 
I u {ui) + I u (u 2 ) 



Ul + u 2 



< Iy (V^iwa) < ^^("1)^(^2) < 



' Ml + U 2 



Moreover, the second and third inequalities hold true for all v > — 1, and t/ie fourth inequality holds true 
for all v > 1/2. In each of the above inequalities equality hold if and only if u\ — u 2 . 

We recall that part (b) of Theorem [T] was proved for v > by Gronwall [27J. Notice also that recently 
Baricz (17) in order to prove the right-hand side of ([4]) proved implicitly part (b) of Theorem [T] For 
reader's convenience we recall below that proof. Moreover, we give a somewhat different proof of this 
part, and two other completely different proofs. 

We note that part (c) of Theorem [1] improves the result of Sun and Baricz [3S] , who proved that 
the function u H> ul v (u) is log-concave on (0,oo) for all v > 1/2. Recently, Baricz and Neuman [TH] 
conjectured that the modified Bessel function I u of the first kind is strictly log-concave on (0, 00) for all 
v > 0. As far as we know, this conjecture is still open and the much sharper result of this kind is of part 
(c) of Theorem [T] 

Proof of Theorem [TJ First we prove the monotonicity and log-concavity properties stated above. 

(a) Recall that the modified Bessel function of the first kind I v is a particular solution of the second- 
order differential equation ([2]) and thus 

(7) I'l{u) = (1 + v 2 /u 2 )I u {u) - (!/«)/». 



Using ([7]) and the left-hand side of ([4]), we obtain that for all u > and v > 1 

_d 

du 



~ul' v (u)~ 




1 




u 2 + v 2 - 2 


'ul' v (u)~ 


2" 


< 


1 




[ i 2 M _ 




ul v (u) _ 




[ I„(u) . 




ul v (u)_ 





1 - V 



< 0. 



(b) Consider the Turanian 

A„(u) = 

which in view of the recurrence relations 

I v -x{u 
and 

I v +\{u) = -{y/u)I v {u) +I' v {u) 



I 2 (u) - I v -x{u)I v+ x{u), 
= {v/u)I v (u) + I' v (u) 
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can be rewritten as follows 



Using ([7]) we get 



A>) = (l + ^ 2 / u 2 )/ 2 ( M )-[/>)] 2 . 
l,o,, \ul' v (u) 



A» = -I 2 (u) 
u 



I V {U) 



It is known (see [IHIHZ]) that the Turan-type inequality A„(u) > holds for all u > and v > — 1, and 
hence the required result follows. We may note incidentally that the result of this part actually follows 
also from the right-hand side of (jj). More precisely, it is easy to see that the function u n- ul' v {u) / I v {u) 
satisfies the differential equation uv'{u) = u 2 + v 2 — v 2 (u), and using the right-hand side of (j4]) it is clearly 
strictly increasing on (0, 00) for all v > — 1. It is important to add here that in fact the right-hand side 
of (@} and the Turan-type inequality A u (u) > are equivalent (see [HI HZ]). 

A third proof of this part can be obtained as follows. By using the infinite series representation of the 
modified Bcsscl function of the first kind we just need to show that the function 



u 1 — ^ 



ul' v {u) 



E 



,2n 



(2n + v){u/2 
n!r(i/ + n+ 1) 



E 

n>0 



(ti/2) 



2/i 



n\T{v 



1) 



is strictly increasing on (0, 00) for all v > —1. To do this let us recall the following well-known result (see 
[201 136) ): Let us consider the power series /(it) = clq + a%u + . . . + a n u n + . . . and g(u) = bo + b\u + . . . + 
b n u n + . . ., where for all n > integer a„£l and b n > 0, and suppose that both converge on (0, 00). If 
the sequence {a n /b n } n >o is strictly increasing, then the function u 1— > f(u)/g(u) is strictly increasing too 
on (0, 00). We note that we can see easily that the above result remains true in the case of even functions. 
Thus, to prove that u M- ul' v (u)/l u (u) is indeed strictly increasing it is enough to show that the sequence 
{ct n }n>o, defined by a n = 2n + v for all n > 0, is strictly increasing, which is certainly true. 
Finally, a fourth proof is as follows. By using the Weierstrassian factorization 

u" T-r / u 2 > 

n>l 



where v > — 1 and j Uin is the nth positive zero of the Bessel function J u of the first kind, we obtain that 
'ul' v {u)' 



_d 

du 



d 
d/t 



2V 



n>l 



Jy.t 



4 E' 

n>l 



2 , ? -2 



> 



for all u > and 1/ > — 1. We note that this proof reveals that the function u H> ul' v {u)/l v {u) is in fact 
strictly decreasing on (— 00, 0) for all v > — 1. This is in the agreement with the fact that the function 
u 1 y ul' v (u) I I v {u) is even, as we can see in the above series representations. 

(c) Owing to Duff |24) it is known that the function u > ^fuK u (u) is strictly completely monotonic, 
and consequently (see [HI p. 167]) strictly log-convex on (0, 00) for each \v\ > 1/2. On the other hand, 
due to Hartman [28] the function u > uI v (u)K u (u) is concave, and consequently log-concave on (0,oo) 
for all v > 1/2. Since u M> 2uI 1 / 2 {u)K 1 / 2 (u) = 1 — e~ 2 " is concave on (0, 00), we conclude that in fact the 
function it t— > uI v {u)K v (u) is concave, and hence log-concave on (0, 00) for all v > 1/2. Now, combining 
these results, in view of the fact that the product of log-concave functions is log-concave, the required 
result follows. 

(d) Using gj and © we obtain that 



_d 

dit 



\u 2 r v {uy 


1 


[ i 2 M _ 





ul' v {u) 

Iu{u) 



- 2 



ul' v (u) 



1 2 



Iu{u) 



< 



+ v 2 + \J u 2 + v 2 - 2 



v + l 

for all u > and v > 0. Observe that the last expression is nonpositive if and only if we have 

2 



v- 1 



2v z 



v + U \ v + l 

A computation shows that this is satisfied if 



- 1 u 2 + v 2 {v 2 - 1) > 0. 
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Now, since v > v o we have 8i/ 3 — 9^ 2 — 2zy — 1 > and thus the proof of part (d) is complete. 

It should be mentioned here that part (a) of this theorem for v > v o actually is an immediate conse- 
quence of this part. More precisely, the proof of part (a) of this theorem can be simplified significantly 
as follows: in view of part (d) of this theorem, the function 

1 u 2 I' v (u) 
u 



u i — y 



is strictly decreasing as a product of two positive and strictly decreasing functions. 

Now, let us focus on the chain of inequalities ©■ To prove this we use Corollary 2.5 from [S]. More 
precisely, the first inequality in ([5]) follows from part (d) of this theorem, while the second inequality 
in (||)]) is an immediate consequence of the fact that I u is a strictly increasing function on (0,oo) for all 
v > — 1. The third inequality in (j6|) means actually the strict geometrical convexity of l v and is equivalent 
to part (b) of this theorem; the fourth inequality is equivalent to part (c) of this theorem. 

Finally, observe that part (a) of this theorem is equivalent to the inequality 

2J„(m)J„(u 2 ) 

Iv(Ui) + I V {U 2 ) 

which holds for all u\,U2 > and v > 1. Moreover, in this inequality equality holds if and only if 

U\ = U2- □ 

The following result is a companion of Theorem [T] for modified Bessel functions of the second kind. 
We note that part (b) of the following theorem is well-known (see for example [25] [3B] I3H] ) , an d part (c) 
was proved by Baricz |17) . For part (b) we give here a different proof, while for part (c) we recall the 
proof from [17] and we present a simple alternative proof. 

Theorem 2. The following assertions are true: 

(a) u i y K' v (u)/K„(u) is strictly decreasing on (0, oo) for all \u\ > 1; 

(b) u i— > K' v (u) I K u (u) is strictly increasing on (0, oo) for all i/£R; 

(c) u i— > uK' v (u) j ' K v {u) is strictly decreasing on (0, oo) for all v 6 M; 

(d) u i— > uK v (u) is strictly increasing on (0, oo) for all v € R; 

(e) u i y u 2 K' u (u) is strictly increasing on (0, oo) for all \v\ > 5/4; 

(f) u n- u 2 K' v (u) is strictly increasing on (2, oo) /or a?Z i/gR, 

In 'particular, for all ui, u 2 > and \i>\ > 1 the following chain of inequalities holds 



(8) 



2K v (u 1 )K v (u 2 ) 
K v (u x ) + K v (u 2 ) 



U\ + u 2 



J < /M^ilM^) < k v (v 7 ^) < 



^(ui) +K v {u 2 ) 



Moreover, the second, third and fourth inequalities hold true for all v G 
U\,U% > £/ie fourth inequality can be improved as 



In addition, for \ v\ > 5/4 and 



(9) 



2miU2 
Ul + U2 



< 



K v (ux) + K v (u 2 ) 



This inequality holds true for all U\,U2 > 2 and v G 
and only if u± = u 2 . 



In each of the above inequalities equality hold if 



Proof. First we prove the monotonicity properties for modified Bessel functions of the second kind. 

(a) Recall that the modified Bessel function of the second kind K v is a particular solution of the 
second-order differential equation @, and this in turn implies that 

(10) K'l{u) = (1 + is 2 /u 2 )K„{u) - (l/u)K' v {u). 

Consequently, by using two times the right-hand side of ([5]), for all u > and v > 1 we have 



d 
du 



\K'M' 




1 




[K 2 (u)_ 




u 


2 K u (u)_ 






< - 




1 











uK{u) 
K v (u) 



uK' v {u) 



uK' v {u) 



uK' v {u) 



K„{u) 



K v {u) 



< 0. 



K v {u) 

On the other hand the function v t-¥ K v (u) is even, and thus from the above result we obtain that indeed 
the function u H> K' v (u) / K 2 {u) is strictly decreasing on (0, oo) for all \v\ > 1. 
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(b) The fact that u h-> K„{u) is log-convex can be verified (see for example [25l |38]) by using the 
Holder- Rogers inequality and the familiar integral representation ([3]), which holds for each u > and 
However, in view of ©, for all n € {0, 1, 2, . . . }, u > and v £ R, we easily have 



Jo 

i.e. the function u t— > K v (u) is strictly completely monotonic. Now, since each strictly completely 
monotonic function is strictly log-convex, we obtain that u i— > K' u {u) / K v {u) is strictly increasing on 
(0, oo) for all 

(c) Consider the Turanian 

Using the recurrence relations 

and 

we have 

K{u) = {l + v 2 /u 2 )Kl{u)-[K' u {u)] 
Combining this with (fTU)) . we obtain [T7] 



A„(u) = Kl(u) - K v ^{u)K v+1 (u). 
tf„_i(u) = -(v/u)K v (u) - K' v {u) 
K y+1 {u) = {u/u)K u (u)-K' u (u) 



K v (u) 

But, the function <—¥ K v (u) is strictly log-convex on R for each fixed u > (see |16j). which implies 
that for all v E R and u > the Turan-type inequality A K (w) < holds. This shows that the function 
u i — ^ uK' v {u)/ 'K v {u) is strictly decreasing on (0, oo) for all ^ G R. Another proof for this part can be 
obtained as follows. First observe that the function u n- uK' v {u) j K v (u) satisfies the differential equation 
uv'{u) = u 2 + v 2 — v 2 (u). On the other hand, it is well-known that K v is strictly decreasing on (0, oo) 
for all v € R. Thus, by using the right-hand side of ([5]) we conclude that u H> uK' v {u) / K v {u) is strictly 
decreasing too on (0, oo) for all v € R. It is important to add here that in fact the right-hand side of ([5]) 
and the Turan-type inequality A„(zt) > are equivalent (see [HI IT7] ). 

(d) By using again the fact that K v is a particular solution of the modified Bessel differential equation, 
i.e. the relation ([1(7)) . we easily have for all u > and v e R 

[uK'Ju)}' = K' v {u) + uK'l(u) = u(l + ^ 2 /"» 2 )^(u) > 0. 

(e) Using (JTOJ) and the left-hand side of ([5]), we obtain 



[u 2 K'M]' = 2 uK(u) u 2 K'l{u) 
K„{u) K„{u) K v {u) 



uK' v {u) 



2 i 2 
U + V 



> u 2 + v 2 - ^Ju 2 vj(y -l) + v 2 



K v {u) 

for all u > and v > 1. The right-hand side of the above inequality is positive if and only if the expression 

Q v {u) = u 4 + [2v 2 - u/{u - l)]u 2 + v 2 {v 2 - 1) 

is positive. It is easy to see that the discriminant of the equation Q v (y/u) — is (5 — Av)v 2 j{v — l) 2 
and this is negative if and only if v > 5/4. Finally, since the function v H > K v (u) is even, the proof is 
complete. 

(f ) In view of Q we obtain that 

/•oo 

u 2 K' v (u) = — w 2 / e~ ucosllt (cosh£)(cosh(i/£))dt 
Jo 

and thus 

/>oo 

[u 2 ^(m)]' =« / (ucosht- 2)e-" cosh *(coshi)(cosh(^))dt > 
Jo 

for all u > 2 and veM. 

Now, let us focus on the inequalities © and (JSJ) . As in the proof of the chain of inequalities ([S]) , we use 
Corollary 2.5 from [5]. The first inequality in © follows from part (a), the second inequality is just the 
strict log-convexity of K v proved in part (b) , while the third inequality is equivalent to the geometrical 
concavity of K v proved in part (c) . The fourth inequality is equivalent to part (d) of this theorem, while 
the inequality (JSJ) is equivalent to part (e). □ 



FUNCTIONAL INEQUALITIES FOR MODIFIED BESSEL FUNCTIONS 



7 



3. Convexity of modified Bessel functions with respect to power means 

In this section we are going to complement and extend the results of the above section. To this aim 
we study the convexity of modified Bessel functions of the first and second kinds with respect to Holder 
means. For reader's convenience we recall here first some basics. 

Let f : [a, b] C R — > R be a strictly monotonic continuous function. The function M v : [a, b) 2 — > [a, b], 
defined by 

+ <p(u 2 ) 



M lp (n 1 ,u 2 ) = p~ 

is called the quasi-arithmetic mean (or Kolmogorov mean) associated to tp, while the function ip is called 
a generating function (or a Kolmogorov- Nagumo function) of the quasi-arithmetic mean M v . A function 
/ : [a, b] C R — > R is said to be convex with respect to the mean M M (or M v — convex) if for all u\,u 2 £ [a, b] 
and all A £ [0, 1] the inequality 

/(AfW(tti,U2)) < M^(f( Ul )J(u 2 )) 

holds, where M^\u\, u 2 ) = (p~ l (\(p(ux) + (1 — X)p(u 2 )) is the weighted version of M v . It can be proved 
easily (see for example 22 ) that / is convex with respect to M v if and only if ip o / o ip~ 1 is convex in 
the usual sense on (f([a, b]). Now, for any two quasi-arithmetic means M v and (with Kolmogorov- 
Nagumo functions <p and ip defined on intervals [a, b] and [c, d\), a function / : [a, b] — ¥ [c,d] is called 
(M v , M^)— convex if it satisfies 

f(MW(u 1: u 2 j) < M^\f( Ul )J(u 2 )) 

for all u\,u 2 £ [a,b] and A £ [0,1], where M^\m,u 2 ) — ip~ 1 (\Tp{ui) + (1 — \)i/j(u 2 )). If the above 
inequality is reversed, then we say that / is (M„, M^)— concave. Due to Aczel [T] it is known from a 
long time ago that if tp is increasing then the function / is (M v , M^)— convex if and only if the function 
ip o f o ip^ 1 is convex in the usual sense on (p([a, b}). This is because, if ip is increasing and we denote with 
s and t the values <p(ui) and ip(u 2 ), then by definition / is [Mm, M^)— convex if and only if 

i/> (/ (tp-'iXs + (1 - X)t))) < A^ (/ (tp-'is))) + (1 - \)ip (/ (t))) 

holds for all s,t £ <p([a,b]) and A £ [0, 1]. See also [22] for more details. 

Now, if ip is decreasing, then clearly the above inequality is reversed, and this in turn implies that the 
function / is {M v , M^)— convex if and only if the function ip o f o p^ 1 is concave in the usual sense on 
</?([a, 6]). Moreover, a similar characterization of (Mm, M^)— concave functions is also valid, depending on 
the monotonicity of the function ip. 

Among the quasi- arithmetic means the Holder means (or power means) are of special interest. They 
are associated to the generating function tp p : (0, oo) — >■ R, defined by 



(p p (u) = 

and have the following form 



u p , ifp^O 
lnu, if p = 0, 



M« K ,) = ( [(l-AK + AnI]V, ifp^O 
y ' ' \ u\u\ A , if p = 0. 

Now, let p and q be two arbitrary real numbers. Using the above definitions of generalized convexities 
we say that a function / : [a, b] C (0, oo) — > (0, oo) is (M Vp , M<p )— convex, or simply (p, q)— convex, if the 
inequality 

(11) /(MW(t»i )U2 )) < MW(/(«i)./(«a)) 

is valid for all p, q £ R, ui, u 2 £ [a, b] and A e [0, 1]. If the above inequality is reversed, then we say that 
the function / is (M m , M v )— concave, or simply (p, q)— concave. Observe that the (1, 1)— convexity is the 
usual convexity, the (1, 0)— convexity is exactly the log-convexity, while the (0, 0)— convexity corresponds 
to the case of the geometrical convexity. We note that motivated by the works [5] [8] and [S], recently Baricz 
[lOj considered the (p,p)— convexity of the zero-balanced Gaussian hypergeometric functions and general 
power series. The (p, q)— convexity of zero-balanced Gaussian hypergeometric functions was considered 
recently by Zhang et al. [33]. 
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The following result gives a characterization of differentiable (p, q)— convex functions and will be applied 
in the sequel in the study of the convexity of modified Bcssel functions of the first and second kinds with 
respect to power means. For a proof see |15) . 

Lemma D. Let p,q 6R and let f : [a, b] C (0, oo) — > (0, oo) be a differentiable function. The function f 
is (strictly) (p,q)— convex ((p,q)— concave) if and only ifut-t u p f' (u)[f (u)] q is (strictly) increasing 
(decreasing) on [a,b]. 

The next result completes and extends parts (a), (b) and (d) of Theorem [1] Notice that if we choose 
in part (b) of Theorem[3]the values p = and q = — 1, then we reobtain part (a) of Theorem[TJ Similarly, 
choosing p = q = in part (a) of Theorem [3] we obtain the strict geometrical convexity stated in part 
(b) of Theorem [TJ Finally, by taking p = q = — 1 in part (b) of Theorem [3] we obtain the monotonicity 
result stated in part (d) of Theorem [TJ 

Theorem 3. Let p, q € K and let v > — 1. Then the following assertions are true: 

(a) if p < and q > 0, then I u is strictly (p,q) — convex on (0,oo); 

(b) if p < and q < 0, then I u is strictly (p, q) — concave on (0, oo) provided if v > —X/q and 

4q(q - l)v 3 - (p 2 - A{q - l))v 2 - 2p 2 v - p 2 > 0; 

(c) if p > and q < —1, then l v is strictly (p, q)— concave on (0, oo) provided if v > 1; 

(d) if p > and q > 0, then L v is strictly (p,q)— convex on (0,oo) provided if v > p/q; 

(e) if p < 1 and q > 1, then I v is strictly (p, q) — convex on (0,oo). 

Proof. For convenience first we introduce the following notation 



_d 

du 



u l - p I' v {u) 



ll- q {u) 



P 



ul' v {u) 



-(1-9) 



ul' u {u) 



We note that in view of Lemma|4]the (p, q)— convexity ((p, q)— concavity) of I v depends only on the sign 
of the expression \ p q v {u). 

(a) This follows easily from the fact that if v > — 1, p < and q > 0, then \ P: q :V (u) > for all u > 0. 
More precisely, from the right-hand side of (U) we have 



yP+l 





~ul' v (u)~ 




'ull(u)~ 


2" 


-p 


[ I v {u) _ 


+ q 


[ I V (U) _ 





> 



for all f > — 1, p < 0, q > and u > 0. It should be mentioned here that this part follows actually from 
part (b) of Theorem [TJ Namely, the function u n- u i-P/^u) [I„(u)} is strictly increasing on (0, oo) 
for all p < 0, q > and v > — 1 as a product of the strictly increasing functions u i-> ul' v (u) / I v (u) and 
u i-> u~ p LQ(u). Now, since for p = q = this part reduces to part (b) of Theorem [TJ the above remark 
reveals that in fact part (b) of Theorem [JJ and part (a) of Theorem [3] are equivalent, 
(b) First assume that p <C and q <C 0. Then by using (|^|) we obtain that 

i q M 



,p+l 



p\fx 



(1 



and this is nonpositive if 



p 2 (u z + ir)< qv l + 



qv 



1 



0<QAu 2 ), 



where Q v (u) = au 2 + bu 



c with v > —1/q, 

2 



qv 



1 



1 , 

This gives a necessary condition to be b 2 
is equivalent to the inequality 

,3 



b = 2qv 



2 av 



1 



2 2/2 2 2\ 

— p , C = V [q V — p ). 



■ 4ac < 0. A computation shows that the condition b 2 — 4ac < 



4q(q - - (p 2 - 4(3 - l))v 2 - 2p 2 v - p 2 > 0. 
Now, assume that p = and q < 0. Then from the left-hand side of ([4]) we have 



v 2 -{l-q) 



Iu{u) 



< 



V v + 1 



qv 



< 
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for all v > — 1/q, q < and u > 0, as we requested. 

(c) This follows directly from part (a) of Thcorcm[TJ More precisely, it is easy to see that the function 
u i y u 1 ~ p I' v {u)I^ l {u) is strictly decreasing on (0, oo) for all v > 1 as a product of the strictly decreasing 
function u h- > ul' v (u)j I%(u) and the decreasing function u n- u~ p IQ +1 (u). Since part (c) of Theorem [3] 
reduces to part (a) of Theorem Q] when p = and q = — 1, the above proof reveals that in fact part (c) 
of Theorem [3] is equivalent to part (a) of Theorem [T] 

(d) Recall that part (b) of Theorem Q] states that I v is strictly geometrically convex on (0, oo) for 
all v > —1, i.e. the function u > ul' v (u) / I„(u) is strictly increasing on (0, oo) for all v > —1. To prove 
that I u is strictly (p, q)— convex on (0, oo) for all p > 0, q > and v > p/q in what follows we show that 
the function u i-» u 1_p i£(u)J3 _1 (u) is strictly increasing as a product of the strictly increasing functions 
u H> ul' u (u) I 'I v (u) and u H> u~ p I%(u). On the other hand, observe that since u i-> ul' v (u) / I v {u) is strictly 
increasing on (0, oo), we obtain that 

ul' v (u) / I v {u) > v 

for all ^ > — 1 and u > (actually for v > this inequality follows directly from the left-hand side of 
Here we used that if u tends to zero then ul' v (u) / I v {u) tends to which can be verified from (@| 
or from 

„,2 



ul' v {u) 

Iu(u) 



The above inequality implies that 



d 






du 


vP 


u p+l 



-p + q 



n>l 



■ 7 2 



> 



l f(-p + ^)>0, 



and with this the proof of this part is complete. 

(e) This follows from the fact that /„ is strictly increasing and convex on (0,oo) for all v > — 1. 
Namely, the function u H> u x ~ p I' v {u)I'i~ 1 (u) is strictly increasing as a product of the strictly increasing 



function u i— > I' v (u) and the increasing functions u i 



□ 



Now, we are going to present the analogous result of Theorem [3] for modified Bessel functions of the 
second kind. We note that part (c) of Theorem [4] (when p = 1 and q = —1) reduces to part (a) of 
Theorem [21 part (e) of Theorem 2] (when p = 1 and q = 0) becomes part (b) of Theorem [21 part (b) of 
Theorem 2] (when p = q = 0) reduces to part (c) of Theorem [51 and part (d) of Theorem 2] (when p = 
and g = l) becomes part (d) of Theorem Finally, observe that if we choose p = — 1 and g = 1 in part 
(a) of Theorem [H then we obtain part (e) of Theorem [H 

Theorem 4. Let p, q G K and Zei i/gl. XTien i/ie following assertions are true: 

(a) j/p < cmd q > 1, then K v is strictly (p,q) — convex on (0, oo) provided if v > 1 and 

4(1 - q)p 2 v 2 + A(q - 2)p 2 v + p 2 {p 2 + 4) < 0; 

(b) if p < and q < 0, then K v is strictly (p, q)— concave on (0, oo); 

(c) if p > and q < 0, then K v is strictly (p,q)— concave on (0,oo) provided if \v\ > —p/q; 

(d) if p > and q > 1, t/ien -RTj, is strictly (p,q)— convex on (0,oo); 

(e) i/p > 1 and q > 0, then K v is strictly (p,q) — convex on (0,oo). 



Proof. For convenience first we introduce the following notation 



l^p,q,v (n) 



_d 

du 



rn 1 " 



Kl- q {u) 



,p+i 



uK(u) 



K v {u) 



(1-9) 



uK(u) 
K v {u) 



Observe that in view of Lemma[3]the (p, q)— convexity ((p, q)— concavity) of K v depends only on the sign 
of the expression /j, Pjgj „(u). 

(a) Notice that for all v e M fixed when u tends to zero uK' v {u) / K v {u) tends to —v. This can be 
verified for example from the integral representation ([3]). On the other hand, in view of part (c) of 
Theorem [5] the function u n- uK' v (u)jK v (u) is strictly decreasing on (0, oo) for all v £ R, and this in 
turn implies that for all v € K and u > the inequality 



(12) 



uK' v (u) / K v (u) < -v 
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holds. We note that actually this follows also from the right-hand side of 
the left-hand side of ([5]) we obtain that 



Now, by using (fT2|) and 



,p+i 



P 



1 



Dv 2 



and the right hand side of the last inequality is nonnegative if and only if 



Q v (u) = u 4 + ( 



2qv 1 - 



+ v 2 (q 2 v 2 -p 2 ) > 0. 



Now, under assumptions the discriminant of the quadratic equation Q v {\/u) — 0, i.e. 

9 

[4(1 - q)p 2 v 2 + A{q - 2)p 2 v + p 2 (p 2 + 4)] 



is negative and with this the proof of this part is complete. 

(b) This follows from the fact that if v € M and p,q < 0, then fi P! q !V (u) < for all u > 0. Namely, 
from the right-hand side of (JSJ we have 



AW,* («) < 



KM 
u p+i 





\uK' v {uY 




\uK' v { U y 


2" 


-p 


[ K u {u) _ 


+ q 


[ K v {u) 





< 



for all v e R, p, q < and u > 0. Here we used that K v is strictly decreasing on (0, oo) for all 
We note here that this part follows actually from part (c) of Theorem [2J Namely, the function u H > 
u l ~ p K' v (u) [K v (u)] q is strictly decreasing on (0, oo) for all p, q < and v g R as a product of the 
strictly decreasing and negative function u n- uK' v (u) / K v (u) and the strictly increasing and positive 
function u u~ p 'K^(u) . Now, since for p = q = this part reduces to part (c) of Theorem [2J the above 
remark shows that in fact part (c) of Theorem [5] is equivalent to part (b) of Theorem 2] 

(c) By using fTS]) and the right-hand side of §5§ we have for all u > 0, p > 0, q < and v > —p/q 

2" 



Hp,q,,A U ) < 



< 



km 

u p+l 

_km 

u p+l 
-(p + qv) 



-P 



uKM 
K v {u) 



uK' v {u) 



uK' v {u) 
K v (u) 

km 

,,p+i 



p-q 
uK' v {u) 



K v {u) 
uK' v (u) 



K v {u) 
< 0. 



K v {u) 

(d) Since p > and q > 1, the function u i-> u~ p K^ -1 {u) is decreasing on (0, oo) for all v e 

9-1 



Now, 

by using part (d) of Theorem [5] we conclude that u i-> v}~ p K' v {u) [K v {u)\ q ~ L is strictly increasing as 
a product of the strictly increasing and negative function u i— > wif ^ (u) and the decreasing and positive 
function u i-> u~ p i4r^ _1 (u). Observe that since for p = and 9=1 this part reduces to part (d) of 
Theorem [21 in fact they are equivalent. Finally, we note that the proof of this part can be obtained also 
simply from the fact that under assumptions Pp, q , u {u) > 0. 

(e) The proof of this part is very similar to the proof of part (d) above. Under assumptions the 
function u H> u 1 ^ v K'j / (u) is decreasing. Consequently, by using part (b) of Theorem [21 the function 
u i v v}~ p K' v (u) [K u (u)] q ~ is strictly increasing as a product of the strictly increasing and negative 

~ p Kl(u). Observe that since 

□ 



function u i— > K' v (u) / K v (u) and the decreasing and positive function u^n 1 
for p = 1 and q — this part reduces to part (b) of Theorem [21 in fact they are equivalent. 

4. Application to the log-concavity of the gamma-gamma distribution 

(0, oo) of the three parameter gamma-gamma 



The probability density function f a ,b,a '■ (0, oo) 
random variable is defined by (see [21]) 



fa,b,a( U ) 



2{ab)- 



2 U 2 



-K a _ 




T{a)T{b)a- 

whcrc a, b > are the distribution shaping parameters, K v stands for the modified Bessel function of 
the second kind, and a > is the mean of the gamma-gamma random variable. The gamma-gamma 
distribution is produced from the product of two independent gamma random variables and has been 
widely used in a variety of applications, for example in modeling various types of land and sea radar 
clutters, in modeling the effects of the combined fading and shadowing phenomena, encountered in 
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the mobile communications channels. Of particular interest is the application of the gamma-gamma 
distribution in optical wireless systems, where transmission of optical signals through the atmosphere is 
involved. For more details see [2TI I23j . 

Now, consider the functions f a ,b.a ■ (0, oo) — > (0,oo) and F a ^, a : (0,oo) — > (0, 1) defined by 



fa,b,a(u) — faM,, 



( cm 2 
\~4ab 



)3-(a+b) 



(ab)u 



a+b-2 



aT(a)T(b) 



-K a _ b (u) 



and 



F a ,b,a{u) 



fa,b.a{t) dt 



T(a)T(b) 



n 2 > 1 

b l,3 



a b 



1 

a, b, 



where G-y 3 is a Meijer G— function [26, ( eq. 9.301]. Here f a .b,a is just a transformation of the proba- 
bility density function f a ,b,a, while F at b, a is the cumulative distribution function of the gamma-gamma 
distribution. 

In probability theory usually the cumulative distribution functions does not have closed-form, and 
thus sometimes it is quite difficult to study their properties directly. In statistics, economics and in- 
dustrial engineering frequently appears some problems which are related to the study of log-concavity 
(log-convexity) of some univariate distributions. An interesting unified exposition of related results on 
the log-concavity and log-convexity of many distributions, including applications in economics, were com- 
municated by Bagnoli and Bergstrom 7 . Some of their main results were reconsidered by Andras and 
Baricz [H] by using the monotone form of l'Hospital's rule. Moreover, by using the idea from [5], recently, 
Baricz |15) showed, among others, that if a probability density function is geometrically concave then the 
corresponding cumulative distribution function will be also geometrically concave. In this section we use 
this result to prove that the cumulative distribution function F a ^, a is strictly log-concave on (0,oo) for 
all a, b, a > 0. This result may be useful in problems of information theory and communications. 

Theorem 5. Let a, b, a > 0. Then the following assertions are true: 

(a) u <— > uf' a b a (u)/ f a ,b,a{u) is strictly decreasing on (0, oo); 

(b) u i y uf' a h a (u) / fa t b,a(u) is strictly decreasing on (0, oo); 

(c) u M> uF' a b a {u) / F at b ta (u) is strictly decreasing on (0, oo); 

(d) a4 F' ab a (u) j -Fa, &,«(«) is strictly decreasing on (0, oo). 

Proof, (a) From part (c) of Theorem [5] we have that the function 

^ U f'a^ a ( U ) __ , j. , U K-b( U ) 



b-2 



K a - b (u) 



fa,b,a{u) 

is strictly decreasing on (0, oo) for all a, 6, a > 0. 

(b) Observe that part (a) of this theorem actually means that the function f a ,b,a is strictly geomet- 
rically concave, i.e. for all a, 6, a > 0, A € (0, 1) and u±, ui > 0, u\ ^ U2 we have 



fa,b,a 

{u\u\ A ) > jf ,i,,a(«l) \ fa,b 



1-A 



Now, changing in the above inequality Ui with 2yJ abui/ a, where i G {1,2}, we obtain 

fa,b,a 

{u^u\ A ) > [fa,b : a{ui)} [fa : b,a(u2)} 

for all a, b, a > 0, A € (0,1) and Ui,it2 > 0, u\ ^ ui- This means that the function f a ,b,a is strictly 
geometrically concave and hence the function u i— ¥ uf b a (u)/ f a ,b,a(u) is strictly decreasing on (0,oo). 

(c) This follows from part (b) of this theorem. Namely, it is known (see [15]) that if the probability 
density function is strictly geometrically concave, then the corresponding cumulative distribution function 
is also strictly geometrically concave. 

(d) Part (c) of this theorem states that the cumulative distribution function F a ^. a is strictly geo- 
metrically concave. Now, by using the fact that f a ,f),a, as a distribution function, is increasing, for all 
a, b, a > 0, A € (0, 1) and u%, U2 > 0, u\ ^ U2 we have 

F a ,b,a (AUI + (1 - X)U 2 ) > F aAa (u*ul- X ) > [F aAa { Ul )] X [Fa^iui)} 1 ^ , 



that is, F at b,a is strictly log-concave on (0,oo) 



□ 
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5. Open Problems 

In this section our aim is to complement the results from the previous sections and to present certain 
open problems, which may be of interest for further research. 

Recall that Neuman [35] proved that the modified Bessel function I v is strictly log-convex on (0, oo) 
for all v G (—1/2,0]. Since I_\i<i{u) — \J ir/ (2m) cosh it, we conclude that in fact /„ is strictly log-convex 
on (0, oo) for all v G [—1/2,0]. Thus, for all v G [—1/2,0] and U\,U2 > the third inequality in ([5) can 
be improved as follows 



Iv {\fU\U 2 ) < Iv ( Ul 2 U2 ) < \/ Iv{ui)I v {U2) 



Moreover, this implies that the function I v is strictly (p,q)— convex on (0, oo) for all v G [—1/2,0], 
p < 1 and q > 0. This can be verified by writing the function u t— > u p I' u (u)I£ _1 (w) as a product of the 
functions u M- I' v (u)/I v {u) and u i-» m 1_p /J(w). 

Concerning Theorem [T] we have the following open problem. 

Question 1. What can we say about the monotonicity of the functions u i— > ul' v (u)/ I 2 (u) and u *— > 
u I' v (u)/I„(u) for \v\ < 1 and v G (— l,fo)j respectively 1 ? Is it true that u i— > \ful v (u) is strictly log- 
concave on (0, oo) for all v > 0? 

Now, concerning Theorem [2j |3j and |4] we may ask the following. 

Question 2. What can we say about the monotonicity of u i— > K' v {u)/K 2 (u) when \v\ < 1? 

Question 3. What can we say about the (p,q)— convexity (concavity) of I v when p > 0, q G (—1,0)? 
Moreover, the conditions for v in parts (b), (c) and (d) of Theorem^ can be relaxed? 

Question 4. What can we say about the (p,q)— convexity (concavity) of K u when p < 1, q £ (0,1)? 
Moreover, the conditions for v in parts (a) and (c) of Theorem^ can be relaxed? 

It is well-known that the function v i-> K v (u) is strictly log-convex on K for all u > fixed (see [T5] ) . 
On the other hand v H> K v (u) is strictly increasing on (0,oo) for all u > fixed. Clearly these imply 
that the function v i— > K v (u) is strictly (p, g)— convex on (0, oo) for all p < 1 and gr > 0, and all fixed 
u > 0. This suggest the following. 

Question 5. What can we say about the (p , q)— convexity (concavity) of the function v t— > K v [u) on 
(0, oo) when p and q are arbitrary real numbers? 

Similarly, the function v H> I v (u) is strictly log-concave on (— l,oo) for all u > fixed (see [H]). On 
the other hand v H> I v (u) is strictly decreasing on (—1, oo) for all u > fixed. Clearly these imply that 
the function v t— > I u {u) is strictly (p, g)— concave on (0,oo) for all p > 1 and g > 0, and all fixed u > 0. 
Thus, it is natural to ask the following. 

Question 6. What can we say about the (p, g) — convexity (concavity) of the function v > Iv(u) on (0, oo) 
when p and q are arbitrary real numbers? And what about the (p, q)— convexity (concavity) of ' v i-> l v (u) 
on (—1, oo) ? 

Due to Laforgia it is known that K' v {u)/K v {u) < —v/u — 1 for all u > and v G (0, 1/2). First 
observe that the above inequality is valid for all v G [0,1/2]. Since Kq(u) = —K\{u) for v = the 
above inequality is equivalent to K\(u) > Kq(u), which is clearly true, since the function v v^r K v (u) 
is strictly increasing on (0,oo) for all u > fixed. Now, since K 1 / 2 (u) — y/ir/ (2u)e~ u we obtain that 
in Laforgia's inequality for v = 1/2 we have equality and since v i— > K u (u) is even, we deduce that 
K' v (u)/ K v {u) < —v/u — 1 holds true for all u > and \v\ < 1/2, with equality for v = 1/2. 

By using this result we obtain that 



[u*K' v (u)]' _ 2 uK' v {u) , u 2 K(u) 



K v (u) K v {u) K v {u) 



UK' (u) 2 2 



< u 2 — u + v 2 — v < 



K v {u) 

for all it G (0, 1) and \v\ < 1/2. This implies that the function u H> u 2 K' v (u) is strictly decreasing on 
(0, 1) for all \v\ < 1/2, i.e. the modified Bessel function of the second kind K v is strictly (—1, 1)— concave 
on (0, 1) for all \v\ < 1/2. This completes parts (e) and (f) of Thcorcm[21 
Taking into account the above discussion we may ask the following. 

Question 7. Is it true that u i— >• u 2 K' v {u) is strictly decreasing on (0,2) for all \u\ < 1/2? 
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In reliability analysis it has been found very useful to classify life distributions (i.e. distributions 
of which cumulative distribution function satisfies F(u) — for u < 0) according to the monotonicity 
properties of the failure rate. By definition a life distribution (with probability density function / and 
survival or reliability function F) has the increasing failure rate (IFR) property if the function u i— >• 
f(u)/F(u) is increasing on (0, oo). Since by definition F(u) = 1 — F(u) for all u > 0, clearly we have 
F (u) = —f(u) for all u > 0. Thus, a life distribution is IFR if and only if u n- — F (u)/F(u) is increasing 
on (0, oo), i.e. the reliability function F is log-concave. It is well-known that if a probability density 
function is log-concave then this implies that the corresponding cumulative distribution function and 
the complementary cumulative distribution function (or survival function) have the same property (for 
more details see [5J [7J [H]). Another class of life distributions is the NBU, which has been shown to be 
fundamental in the study of replacement policies. By definition a life distribution satisfies the new-is- 
better-than-used (NBU) property if u H> log F(u) is sub-additive, i.e. 

F(ui + u 2 ) < F( Ul )F(u 2 ) 

for all u%,U2 > 0. The corresponding concept of a new-is-worse-than-used (NWU) distribution is defined 
by reversing the above inequality. The NBU property may be interpreted as stating that the chance 
F(m) that a new unit will survive to age u\ is greater than the chance F(ui+u 2 )/F(u 2 ) that an unfailed 
unit of age u 2 will survive an additional time u\. It can be shown easily that if a life distribution is IFR 
then it is NBU (see for example [H]), but the inverse implication in general does not hold. Since the 
most important life distribution satisfies the NBU property it is natural to ask the following. 

Question 8. Is it true that the gamma-gamma distribution satisfies the NBU property? 

To answer this question it would be enough to prove that the probability density function f a .b,a is log- 
concave, and for this in view of part (b) of Theorem[S]it is quite enough to show that f a ,b,a is increasing. 
Similarly, observe that for the log-concavity of f a ,b,a we just need to show that f a ,b,a is increasing and 
log-concave. However, by part (a) of Theorem [5] if f a ,b,a is increasing, then it is log-concave. Thus, 
to prove that the gamma-gamma distribution is NBU we need to show that either f a ,b,a or f a ,b,a is 
increasing. 
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